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Using  accurate  band  structures  of  InAs,  InSb,  and  two  Hg^CcfTe  alloys,  we  calculate  the  change  in  refractive 
index  caused  by  the  photoexcited  electrons  and  holes.  The  effects  of  both  free-carrier  absorption  (FCA)  and 
one-photon  absorption  are  considered.  We  find  that  the  change  in  refractive  index  varies  nonlinearly  with  the 
density  of  photoexcited  carriers  and  that  the  generally  neglected  FCA  contribution  is  significant  in  InAs,  owing 
to  its  weak  spin-orbit  coupling.  ©  2006  Optical  Society  of  America 
OCIS  codes:  190.5970,  160.4760. 


1.  INTRODUCTION 

Semiconductor  materials  such  as  InSb,  InAs,  and 
HgCdTe  are  used  in  optoelectronic  applications  including 
quantum-cascade  lasers,  terahertz  generation,  and  satu¬ 
rable  absorbers.1,2  These  applications  involve  generation 
of  high  levels  of  charge  carriers  that  render  it  important 
to  understand  the  effect  of  photoexcited  carriers  on  the 
material  optical  properties.  In  addition,  high-intensity 
light  propagation  in  semiconductors  can  be  strongly  af¬ 
fected  by  the  generation  of  a  high  density  of  photoexcited 
carriers.  The  incident  light  creates  electron-hole  pairs, 
some  of  which  recombine  through  a  variety  of  mecha¬ 
nisms  including  the  Auger,  radiative,  and  Shockley- 
Read-Hall  processes.  Not  all  electron-hole  pairs  will  re¬ 
combine  before  the  entry  of  the  next  stream  of  photons, 
and  the  excess  excited  carriers  will  alter  the  photon  ab¬ 
sorption  both  above  and  below  the  fundamental  gap.  As  a 
result,  the  refractive  index  changes  from  its  value  at  equi¬ 
librium.  Although  there  are  a  number  of  studies  to  accu¬ 
rately  calculate  the  absorption  and  carrier  relaxation  in 
semiconductors,  only  simplified  models  that  use  effective 
mass  band  structures,  the  Boltzmann  distribution  for  the 
carriers,  and  one-photon  absorption  (OPA)  across  the 
bandgap  are  currently  being  employed  in  the  evaluation 
of  refractive  index  changes.3-8  These  calculations  predict 
that  the  change  in  refractive  index  varies  linearly  with 
the  carrier  density  and  increases  with  lattice  tempera¬ 
ture.  The  contribution  from  free-carrier  absorption  (FCA) 
is  considered  only  within  the  Drude  approximation.9 

In  this  paper,  we  use  full  band  structures  and  Fermi- 
Dirac  (FD)  statistics  for  photoexcited  carriers  to  system¬ 
atically  calculate  the  change  in  refractive  index  arising 
from  both  OPA  and  FCA  in  four  direct  bandgap  semicon¬ 
ductors:  two  having  bandgaps  less  than  0.21  eV  (InSb  and 
Hgo.762Cdo.238Te)  and  two  having  a  bandgap  of  approxi¬ 
mately  0.35  eV  (InAs  and  Hg0  656Cd0.344Te).  The  results 
from  this  accurate  calculation  differ  from  the  existing  re- 
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suits  not  only  in  magnitude,  but  also  in  the  trend.  The  cal¬ 
culated  change  in  retractive  index  varies  highly  nonlin- 
early  with  the  excess  carrier  density  in  the  two  smaller 
bandgap  materials  but  varies  almost  linearly  in  the  two 
larger  bandgap  materials  and  has  a  considerable  contri¬ 
bution  from  FCA  in  InAs,  where  the  spin-orbit  (SO)  split¬ 
ting  is  small. 

2.  FORMALISM 

The  photoexcited  carriers  in  the  conduction  and  valence 
bands  modify  the  absorption  spectrum,  which  results  in  a 
change  in  the  refractive  index  through  the  Kramers- 
Kronig  relation  between  the  real  part  of  the  refractive  in¬ 
dex  n i  and  the  absorption  coefficient  a(co): 

c  f00  a(co')dco' 

n1(co)-l  =  -P\  —2 - 2",  (D 

77  J0  CO2 -CO2 

where  c  is  the  speed  of  light,  co  is  the  angular  frequency, 
and  P  indicates  the  principle  value  of  the  integral.  The  re¬ 
fractive  index  change  A ri\  arising  from  the  excess  carriers 
can  be  calculated  from10 

c  f00  a(co' ,AN)  -  a(co' ,0) 

A  n^co)  =  -P  - ^ - 2 - do/  ’  (2) 

77  J0  CO2 -CO2 

where  a(co',AN)  and  a(co' ,  0)  are  the  absorption  coeffi¬ 
cients  at  angular  frequency  of  co'  in  the  presence  of  pho¬ 
toexcited  electron-hole  pairs  with  the  density  of  AN  and 
in  equilibrium,  respectively. 

The  absorption  coefficient  is  related  to  the  imaginary 
part  of  the  dielectric  function,  e2(co),  via  a(co) 
=  [a)62(oj)]/[^1(w)c].  For  a  small  change  in  refractive  index, 
n^co)  in  this  relation  can  be  replaced  by  n°(co),  which  is 
the  refractive  index  for  the  system  in  equilibrium.  The  ab¬ 
sorption  change  can  then  be  written  in  terms  of  the 
change  in  the  dielectric  function, 
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0) 

a(u,AN)  -  =  o  [e2(to,AiV)  -  e2(w,0)],  (3) 

n{(co)c 

where  e2(o),0)  and  e2(a),AN)  are  the  imaginary  part  of  the 
dielectric  function  for  the  system  in  equilibrium  and  in 
the  presence  of  additional  electron-hole  pairs  with  den¬ 
sity  AN,  respectively. 

Among  many  absorption  processes,  the  FCA  and  OPA 
are  most  sensitive  to  the  carrier  occupation  in  a  small 
bandgap  material.  The  OPA  excites  an  electron  from  the 
valence  band  into  the  conduction  band  when  the  photon 
has  energy  greater  than  the  bandgap.  To  achieve  this  in¬ 
terband  transition,  the  initial  state  in  the  valence  band 
must  be  occupied  and  the  final  state  in  the  conduction 
band  must  be  vacant.  The  photoexcited  carriers  change 
the  occupation  in  both  conduction  and  valence  bands 
when  compared  with  the  system  in  equilibrium  and  thus 
modify  the  OPA.  When  the  energy  of  a  photon  is  not  large 
enough  to  enable  valence-band  to  conduction-band  transi¬ 
tions,  the  carriers  can  still  absorb  the  photon  through  in¬ 
traband  transitions  (with  the  help  of  phonons)  and 
intervalence-band  transitions,  commonly  known  as  FCA. 
Since  both  these  processes  are  proportional  to  the  number 
of  carriers,  the  photoexcited  carriers  will  modify  the  FCA 
spectrum  as  well. 

We  calculate  the  change  in  the  dielectric  function  aris¬ 
ing  from  the  OPA  and  FCA  by  electrons  and  holes  using 
accurate  band  structures  and  wave  functions  from  a  long- 
range  tight-binding  Hamiltonian.11  The  OPA  contribution 
to  €2(o),AN)  can  be  calculated  from 


1  _  1 

Rnk  =  2j  ^LO  +  ^  +  2 


n  k  ±  q 


72k 


Q 


72k  +  q 


HeiP 


72k  +  q 


(7) 


where  n  is  the  band  index,  ho)LO  is  the  LO  phonon  energy, 
and  7VLc>  is  the  number  of  LO  phonons  that  can  be  ob¬ 
tained  from  the  Bose-Einstein  distribution  function.  We 
have  used  the  Frohlich  Hamiltonian  to  describe  the  cou¬ 
pling  between  the  electrons  and  the  LO  phonons.  The  sec¬ 
ond  part  arises  from  photon  absorption  by  holes,  without 
a  phonon,  making  a  direct  transition  between  the  heavy- 
hole  and  light-hole  valence  bands.  This  contribution  is 
calculated  from 


e2(o),AN) 


4ir2e2 

3 772  q  ft)2 


IPrara' 

mm'  k 


(k)|2/*(Emk)[l-/*(Em,k)] 


X£(£m,k-£mk-  ha>), 


(8) 


where  m  and  m'  are  valence-band  indices. 

It  should  be  noted  from  Eqs.  (2)  and  (3)  that  the  dielec¬ 
tric  function  needs  to  be  computed  over  a  very  large  fre¬ 
quency  range  even  if  At2X  is  evaluated  only  near  the  band- 
edge  frequencies.  In  addition,  a  careful  numerical 
evaluation  of  the  integral  in  Eq.  (2)  is  required  because  of 
the  singularity  in  the  integrand. 


4772C2 

e2(W,AA0  =  — 2  |p„m(k)|2f(£mk)[l  -f(Enk)] 

3  772  o  ft)  ^mk 

x  -  Emk  -  h  o)),  (4) 

where  e  is  the  electron  charge;  7720  is  the  free-electron 
mass;  n  and  m  are  the  conduction-  and  valence-band  in¬ 
dices,  respectively;  En^  and  Emk  are  the  corresponding 
band  energies;  and  pnm  are  the  interband  momentum  ma¬ 
trix  elements,  f  is  the  charge-carrier  distribution  function, 


f(Emk)  =  [  l  + 


f(Enk)  =  [  1+eWam)]-1,  (5) 

where  /2e  and  /jl ^  are  the  Fermi  levels  for  electrons  and 
holes,  respectively,  which  are  determined  from  the  den¬ 
sity  of  photoexcited  electron-hole  pairs,  AN.  When  AN 
=  0,juLe  and  jmh,  are  equal. 

The  FCA  contribution  to  e2(o),AN)  includes  two  parts. 
The  first  is  that  electrons  in  the  conduction  bands  (or 
holes  in  the  valence  bands)  absorb  a  longitudinal-optical 
photon  and  make  a  transition  to  a  higher-energy  state 
within  the  same  band.  This  contribution  is  calculated 
from 


€2((x),AN) 


87r3e4ojLo  /  1 
huPrn^V  \  /Coo 


-ERAk, 

*0/  n~k 


(6) 


3.  RESULTS  AND  DISCUSSION 

To  test  the  accuracy  of  the  band  structures  and  dipole  ma¬ 
trix  elements  used  in  our  calculations,  we  have  computed 
the  OPA  coefficients  in  InSb  and  InAs  and  found  them  to 
be  in  excellent  agreement  with  the  measured  values.  The 
dielectric  function  and  the  change  in  the  refractive  index 
are  then  calculated  for  the  two  materials  InAs  and 
Hgo.656Cd0  344Te  with  bandgap  Eg  of  0.35  eV  and  the  two 
smaller  bandgap  materials  in  Sb  at  E^=0.173eV  and 
Hgo.762^d0  238Te  with  E^  =  0.206eV.  The  spectra  of  n®(a)) 
used  in  our  calculations  are  from  the  literature.12-15 

Figure  1  shows  the  calculated  change  in  refractive  in¬ 
dex  as  a  function  of  wavelength  due  to  the  FCA  and  OPA 
in  InSb  at  two  densities  of  photoexcited  electron-hole 
pairs — AN=  1016  and  1017  cm-3.  We  see  that  the  FCA  con¬ 
tribution  depends  weakly  on  the  wavelength.  However, 
the  OPA  contribution  shows  a  steep  change  near  the 
bandgap.  This  variation  can  be  easily  understood  from 
Eq.  (2).  In  the  presence  of  photoexcited  carriers,  the  band- 
edge  states  are  filled,  and  the  joint  density  of  states  of 
electrons  and  holes  available  for  absorption  decreases  and 
the  OPA  is  reduced.  Thus  the  numerator  of  the  integrand 
in  Eq.  (2)  is  either  zero  (for  o)f  <  o)g)  or  negative  (for  o/ 
>  o)g)  with  o)g  being  the  bandgap  frequency.  When  the  fre¬ 
quency  o)  crosses  cog,  the  denominator  of  the  integrand 
changes  sign  and  consequently  An^oo)  changes  sign. 

The  calculated  value  of  At2X  at  \  =  9.5yL<m  is  shown  in 
Fig.  2  as  a  function  of  photoexcited  carrier  density  in  InSb 
and  Hg0  762Cd0  238Te.  We  see  that  in  both  materials  the 
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Wavelength  [pm] 

Fig.  1.  Refractive  index  change  due  to  the  FCA  and  OPA  as  a 
function  of  the  photon  wavelength  in  InSb  with  AN=1016  and 
1017  cm-3  at  room  temperature. 


AN  [cm'3] 


Fig.  2.  Refractive  index  change  at  X  =  9.5  pm  due  to  the  FCA  and 
OPA  as  a  function  of  A N  in  InSb  and  Hg0  762Cd0  238Te  (HCT)  at 
room  temperature. 

OPA  contribution  dominates  over  the  FCA  contribution  to 
the  A n\.  Interestingly,  the  FCA  contribution  is  positive, 
which  is  opposite  of  the  OPA  contribution  at  this  wave¬ 
length.  Similarly  calculated  results  at  X  =  4.8pm  for  the 
larger  bandgap  materials — InAs  and 

Hgo.656Cd0 .344Te — are  plotted  in  Fig.  3.  In  this  case,  both 
FCA  and  OPA  reduce  the  refractive  index.  Although  both 
materials  have  the  same  bandgap,  in  the  HgCdTe  alloy 
the  dominant  contribution  is  from  the  OPA,  whereas  in 
InAs  both  FCA  and  OPA  have  comparable  contributions. 
To  understand  the  different  levels  of  contribution  of  FCA 
to  the  refractive  index  change,  we  plot  the  FCA  cross  sec¬ 
tion  (rasa  function  of  wavelength  in  Fig.  4  for  these  two 
materials.  The  numerator  of  the  integrand  in  Eq.  (2)  for 
the  FCA  process  is  proportional  to  this  cross  section,  and 
the  FCA  contribution  to  Apx  at  a  photon  wavelength  X  is 
roughly  proportional  to  the  difference  in  the  areas  below 


Yu  et  al. 

and  above  that  X  in  Fig.  4.  We  see  that  the  absorption  for 
wavelengths  below  X  =  4.8  pm  is  negligible  in  InAs,  but  is 
significant  in  the  HgCdTe  alloy.  Owing  to  a  larger  cancel¬ 
lation,  the  FCA  contribution  is  smaller  in  the  HgCdTe  al¬ 
loy.  We  explore  further  to  understand  the  origin  of  the  dif¬ 
ference  in  the  absorption  spectrum  of  Fig.  4.  In  our 
calculations  we  found  that  most  of  the  FCA  arises  from 
the  transitions  of  holes  from  the  heavy-hole  to  light-hole 
bands.  Hence  the  different  FCA  spectra  can  be  under¬ 
stood  from  the  valance-band  structures  in  the  two  mate¬ 
rials.  We  plot  in  Fig.  5  the  energy  difference  between  the 
heavy-hole  and  the  light-hole  valence  bands  as  a  function 
of  wave  vector  k  in  symmetry  directions  for  InAs  and 
Hgo.656Cdo.344Te.  In  both  cases  the  energy  difference  in- 


AN  [cm’3] 

Fig.  3.  Refractive  index  change  at  X  =  4.8  pm  due  to  the  FCA  and 
OPA  as  a  function  of  A N  in  InAs  and  Hg0  656Cd0  344Te  (HCT)  at 
room  temperature. 


Fig.  4.  Cross  section  of  the  FCA  a  as  a  function  of  the  photon 
wavelength  in  InAs  and  Hg0.656Cd0.344Te  (HCT)  at  room 
temperature. 
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Fig.  5.  Energy  difference  between  the  heavy-hole  and  light-hole 
valence  bands  as  a  function  of  wave  vector  for  InAs  and 
Hgo.656Cd0_344Te  (HCT)  at  room  temperature. 


Fig.  6.  Total  change  in  the  refractive  index  A nl  at  A  =  4.8  jmm  for 
InSb  and  Hg0  762Cd0  238Te  (HCT-238)  and  at  A  =  9.5  jmm  for  InAs 
and  Hg0  656Cd0  344Te  (HCT-344)  as  a  function  of  AN  at  room 
temperature. 


Table  1.  Fit  Parameters  to  Obtain  the  Calculated 
Change  in  Refractive  Index  at  300  K 


Parameter 

InAs 

Hgi-xCdjTe 

{x = 0.344) 

InSb 

Hgi-^Cd.Te 
(x= 0.238) 

a 

0.13067 

0.14536 

0.077854 

0.12091 

b 

8.1952 

7.079 

2.9875 

1.6804 

x  10-19 

X  10"19 

x  10-18 

X  10"18 

creases  rapidly  at  small  values  of  k,  but  saturates  at 
0.25  eV  as  the  wave  vector  increases  further  in  InAs.  This 
small  saturation  value  can  be  attributed  to  the  weak  SO 
coupling  (0.38  eV)  in  InAs,  which  repels  the  light-hole  va¬ 
lence  band  closer  to  the  heavy-hole  valence  band.  Absorp- 
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tion  of  photon  with  a  larger  energy  than  this  saturation 
value  is  not  possible.  In  the  HgCdTe  alloy,  however,  a 
larger  SO  coupling  (0.95  eV)  allows  the  light-hole  band  to 
move  farther  away  from  the  heavy-hole  band,  resulting  in 
absorption  at  shorter  wavelengths.  Hence  the  relative 
contribution  from  OPA  and  FCA  depend  critically  on  the 
subtle  details  of  the  band  structures. 

The  sum  of  the  FCA-  and  OPA-induced  change  in  re¬ 
fractive  index,  calculated  at  4.8  jmm  for  the  large  gap  ma¬ 
terials  and  at  9.5  jmm  for  the  small  gap  materials,  is 
shown  in  Fig.  6  as  a  function  of  the  density  of  photoex- 
cited  electron-hole  pairs.  Clearly  the  variation  is  nonlin¬ 
ear  for  the  two  smaller  bandgap  materials.  An  accurate  fit 
to  the  total  change  indicated  that  variation  is  nonlinear 
even  for  the  larger  bandgap  materials.  Excellent  fits  to  all 
four  materials  are  obtained  by  using  a  functional  form 
A^1=a(l-c_6AiV).  Only  for  small  AiV(<1016  cm-3)  does  the 
refractive  index  change  linearly  with  AN.  The  fitted  val¬ 
ues  of  a  and  b  are  given  in  Table  1. 

It  is  interesting  to  estimate  the  effect  of  temperature  on 
the  OPA-induced  change  in  the  refractive  index.  In  our 
calculation  with  FD  statistics  and  full  band  structures, 
the  excess  carrier-induced  change  in  OPA  is  nearly  tem¬ 
perature  independent.  However,  for  example,  the  band- 
gap  in  InSb  changes  from  0.173  eV  at  300  K  to  0.234  eV  at 
77  K.  Since  the  contribution  arising  from  OPA  is  nonzero 
only  when  co'  >  o)g,  the  denominator  of  the  integrand  in 
Eq.  (2)  is  larger  at  T=  77  K  than  at  300  K,  resulting  in  a 
smaller  change  in  the  index.  This  conclusion  is  in  contra¬ 
diction  with  previous  studies,5’6  which  used  simplified 
band  structures  and  a  Boltzmann  distribution. 


4.  CONCLUSIONS 

We  have  calculated  the  change  in  refractive  index  arising 
from  both  the  OPA  and  the  FCA  using  full  band  struc¬ 
tures  for  four  materials:  InSb,  Hg0  762Cdo.238Te,  InAs,  and 
Hgo.656Cd0.344Te.  We  focus  on  the  refractive  index  change 
at  A  =  9.5  jmm  for  the  first  two  materials  that  have  band- 
gaps  less  than  0.21  eV  and  at  A  =  4.8  jmm  for  the  last  two 
materials  that  have  a  bandgap  of  0.35  eV.  The  calcula¬ 
tions  show  that  the  FCA  contribution  is  not  negligible  and 
is  comparable  to  the  OPA  contribution  in  InAs  that  has  a 
small  SO  splitting.  The  calculated  dependence  of  the 
change  in  refractive  index  on  the  density  of  photoexcited 
electron-hole  pairs  is  strongly  nonlinear  in  InSb  and 
Hgo.762Cd0.238Te  and  moderately  nonlinear  in  InAs  and 
Hgo.656Cd0.344Te.  Furthermore,  the  change  in  refractive 
index  for  InSb  due  to  the  OPA  is  smaller  at  lower  tem¬ 
peratures  than  at  room  temperature  because  of  a  larger 
bandgap  at  low  temperatures. 
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